Mori. Not. R. Astron. Soc. 000. ITH2T1 (20021 Printed 25 February 2013 (MN LM^jX style file v2.2) 
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ABSTRACT 

Broad absorption lines (BALs) in quasar spectra indicate high- velocity outflows that 
may be present in all quasars and could be an important contributor to feedback to 
their host galaxies. Variability studies of BALs help illuminate the structure, evolution, 
and basic physical properties of the outflows. Here we present further results from an 
ongoing BAL monitoring campaign of a sample of 24 luminous quasars at redshifts 
1.2 < Z < 2.9. We directly compare the variabilities in the C iv A1549 and Siiv A1400 
absorption to try to ascertain the cause(s) of the variability. We find that Siiv BALs 
are more likely to vary than Civ BALs. When looking at flow speeds >— 20 000 km 
s _1 , 47 per cent of quasars exhibited Si IV variability while 31 per cent exhibited 
C IV variability. Furthermore, ^50 per cent of the variable Si iv regions did not have 
corresponding C iv variability at the same velocities, while nearly all occurrences of 
C iv variability had corresponding changes in Si iv. We do not find any correlation 
between the absolute change in strength in C iv and in Siiv, but the fractional change 
in strength tends to be greater in Siiv than in Civ. When both Civ and Siiv varied, 
those changes always occurred in the same sense (either getting weaker or stronger). 
We also include our full data set so far in this paper, which includes up to 10 epochs of 
data per quasar. The multi-epoch data show that the BAL changes were not generally 
monotonic across the full ~5 to ~8 yr time span of our observations, suggesting that 
the characteristic time-scale for significant line variations, and (perhaps) for structural 
changes in the outflows, is less than a few years. Coordinated variabilities between 
absorption regions at different velocities in individual quasars seems to favor changing 
ionization of the outflowing gas as the cause of the observed BAL variability. However, 
variability in limited portions of broad troughs fits naturally in a scenario where 
movements of individual clouds, or substructures in the flow, across our lines-of-sight 
cause the absorption to vary. The actual situation may be a complex mixture of 
changing ionization and cloud movements. Further discussion of the implications of 
variability, e.g., in terms of the size and location of the outflowing gas, will be presented 
in a forthcoming paper. 

Key words: galaxies: active - quasars: general - quasars:absorption lines. 



1 INTRODUCTION 

High- velocity wind-like outflows are an important part of the 
quasar system and a potential contributor to feedback to the 
host galaxy. These outflows may play a crucial role in the 
accretion process in quasars and the growth of supermassive 
black holes (SMBHs), by allowing the accreting material to 
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release angular momentum. Furthermore, quasar outflows 
may provide enough kinetic energy feedback to affect star 
formation in the quasar host galaxies, to aid in 'unveiling' 
dust-enshrouded quasars, and to hel p distribute metal-rich 
gas to the interg alactic medium fe.g.. |Pi Matteo et alj|2005l : 
iMollet alj|2007f ). 

Broad absorption lines (BALs) are the most prominent 
signatures of accretion disk outflows seen in quasar spec- 
tra. BALs are defined as absorption troughs with velocity 
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_1 at depths >10% below the contin- 
uum IjWevmann et al.lll991f). and they appear in the spectra 



widths >2000 km s 

of ~ 10-15% of quasars (|Reichard et akll2003l ; I Trump et al.l 



2006; Gibson et al.ll2009l ). Since BALs are seen in just a frac- 
tion of quasar spectra, their presence could represent a phase 
in the evolution of a quasar and/or particular orientations 
where the outflow lies between us and the quasar emission 
sources. 

The location and three-dimensional structure of quasar 
outflows are poorly understood. Sophisticated models pre- 
dict these outflows as arising from a rotating accre- 
tion disk, with acceleration to high speeds by radia- 
tive and/or magneto-centrifugal forces (IMurrav et al.lll995l ; 
IProga fc Kallmanl |2004| ; IProgal |2007| ; lEverettl 12005^ . Im- 
proved observational constraints are necessary to test these 
models and to estimate mass-loss rates, kinetic energy yields, 
and the role of quasar outflows in feedback to the surround- 
ing environment. 

One way to obtain constraints on quasar outflows is 
to study the variability in their absorption lines, which 
can provide information on the structure and dynamics 
of the outflowing gas. Two possible causes of this ob- 
served variability are movement of gas across our line of 
sight to the quasar and changes in ionization ([B arlow 
19931 ; I Wise et all 12004 iMjsawa et al.ll2005l; ILundgren et al 



20071 ; JCibson et alj|2008l : lHamann et alj|2008h . Variabilis 
on shorter time-scales can place constraints on the dis- 
tance of the absorbing material from the central SMBH. 
Shorter variability time-scales indicate smaller distances, 
based on nominally shorter crossing times for moving 
clouds (|Hamann et al.l 120081 ; ICapellupo et al.l 120111 ) or the 
higher densities requ ired for shorter recombination times 
IjHamann et al.ll 19971 ). Measurements of variability on longer 
(multi-year) time-scales provide insight into the homogene- 
ity and stability of the outflowing gas. If no variability is 
detected on long time-scales, then this indicates a smooth 
flow with a persistent structure. Overall, results of variabil- 
ity studies provide information on the size, kinematics, and 
internal makeup of sub-structures within the outflows. Fur- 
thermore, variability studies can address the evolution of 
these outflows as absor ption lines have been observed to 
appear and disappear (Hamann et al. 2008; Lcighlv et al. 
120091 ; iKrongold et al-lkoioHRodriguez Hidalgo et aLlbOlU T 
or they can evolve from one type of outflow feature to 
another (e.g., from a mini-BAL to a BAL or vice versa; 
iGibson et al.ll2010l ; Rodriguez Hidalgo et al. in preparation; 
also this work). 

Most of the existing work on BAL variability has fo- 
cused on var i ability in C IV A1549 over two epoch s (e.g., 



Barlowl Il993l; ILundgren et al.l 120071 ; IGibson et all [2008). 
Gibs on et al.ll2008l detected Civ BAL variability in 12 out 
of 13 BAL quasars (92 per cent) over multi-year time-scales. 
None of thes e studies foun d clear evidence for acceleration 
in the BALs. IGibson et al.l (|2010l ) reports on variability on 
multi-month to multi-year rest-frame time-scales, using 3-4 
epochs of data for 9 BALQSOs and found that BALs gener- 
ally do not vary monotonically over time. Their study also 
makes comparisons between variability in Siiv A 1400 ab- 
sorption and variability in C IV, and their results include a 
correlation between fractional change in EW in Si IV and 
Civ. 

This work is the second paper in a series on BAL vari- 



ability. The first paper, ICapellupo et al.l (1201 ll ; hereafter, 
Paper 1), introduced our ongoing monitoring programme of 
a sample of 24 BA L quasa r s. We began with a sample of 
BAL quasars from I Barlowl ((1993), which includes spectra 
of the C IV absorption region and, in most cases, coverage 
of the Si IV absorption region as well. We have re-observed 
these quasars to provide a longer time baseline over which to 
study variability, as well as to obtain multiple epochs of data 
per object. We currently have up to 10 epochs of data per 
quasar up to March 2009, covering rest-frame time intervals 
(At) from 15 days to 8.2 yfl 

Paper 1 focused on a subset of the data from this mon- 
itoring programme to look for basic trends in the data be- 
tween variability and other properties of the absorbers, as 
well as to directly compare short-term and long-term vari- 
ability within the same sample of quasars. Paper 1 took a 
novel approach to studying BAL variability by introducing 
a measurement of BAL strength within portions of a trough, 
instead of using equivalent width (EW) measurements. Pa- 
per 1 discusses variability in just two different time inter- 
vals: a short-term interval of 0.35—0.75 yr and a long-term 
interval of 3.8—7.7 yr. We found that 39 per cent (7/18) of 
the quasars varied in the short-term, whereas 65 per cent 
(15/23) varied in the long-term data. The variability most 
often occurred in just por tions of a BA L trou gh, which is 
similar to the findings of IGibson et al.l (|2008l ). We found 
that the incidence of variability was greater at higher ve- 
lo cities and in weaker p ortions of BAL troughs. Similarly, 
in ILundgren et al.l (|2007l ). the strongest occurrences of BAL 
variability occurred at velocities < — 12 000 km s _1 and in 
features with smaller equivalent widths. Overall, the results 
of Paper 1 are broadly consistent with previous work on 
BAL variability. 

In this paper, we extend the analysis of Paper 1 by look- 
ing at variability in Si IV and comparing it to the variability 
results for C IV. Expanding our study to include Si IV ab- 
sorption can help constrain theories on the cause(s) of BAL 
variability. C and Si have different abundances, if solar abun- 
dances are assumed, and they have di fferent ionization prop- 
erties (e.g. lHamann et ai1l2008l . l201ll ). By examining if Civ 
and Si IV have different variability properties, and how they 
differ, coupled with these differences in abundances and ion- 
ization properties, we can gain new insight into the cause(s) 
of BAL variability. 

Our dataset is uniquely suited to this study because we 
have coverage of the Si IV line for nearly our entire sample 
(22 out of 24 quasars). In addition to the larger sample size, 
we go beyond existing work by adopting a method of mea- 
suring the absorption strength in portions of BALs, instead 
of EW. Equivalent width measurements apply to an entire 
feature and are less sensitive to changes in small portions of 
troughs. Our method of measuring portions of troughs also 
allows more direct comparisons between the behavior of C IV 
and Si IV variability. 

We also include the entire dataset so far to look at vari- 
ability in C IV and Si IV over multiple epochs. This work 
contains up to 10 epochs of data per quasar, and includ- 
ing all of these epochs will provide better insight into the 
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characteristics of BAL variability. Increasing the number of 
epochs provides new information on whether BALs change 
monotonically over time or whether they can vary and then 
return to an earlier state. In Paper 1, we reported that typ- 
ically only portions of BALs varied. Multi-epoch data can 
tell us whether variability only occurs in those specific ve- 
locity intervals or if the velocity range over which variability 
occurs can change over time. We also highlight several indi- 
vidual interesting cases of variability that can further help 
us understand BAL outflows. Section 2 below reviews the 
quasar sample and analysis introduced in Paper 1, Section 3 
describes our results, Section 4 summarizes the results so far 
from Paper 1 and the current work, and Section 5 discusses 
the results and their implications. 



2 DATA AND ANALYSIS 

2.1 Observations and quasar sample 

In this work, we use the same sample of 24 BAL quasars 
introduced in Paper 1. This sample is based on the set of 
BALQSOs studied in iBarlowl (| 19931 ). The sample selection 
and general characteristics are described in Paper 1. These 
data were obtained from the Lick Observatory 3-m Shane 
Telescope, using the Kast spectrograph. Most of the spectra 
we use from that data set have a resolution of R = A/AA « 
1300 (230 km s _1 ). For epochs where this resolution is not 
available, we use data taken at R ~ 600 (530 km s _1 ). BALs 
are defined to have a width of at least 2000 km s - , so 
either of these resolutions is sufficient to measure the lines 
and study their variabilities. The wavelength coverage of 
each spectrum covers at least the Si IV through C IV emission 
lines, and most cover at least the Lya A1216 through C IV 
emission lines. 

We have been re-observing 23 of the BALQSOs from 
IBarlowl (|l993T ) at the MDM Observatory 2.4-m Hiltner tele- 
scope, using the CCDS spectrograph with a resolution of 
R ~ 1200 (250 km s _1 ). The observations used in this work 
were taken in January and February 2007; January, April, 
and May 2008; and January and March 2009. We used the 
same spectrograph setup each time, varying only the wave- 
length range in order to observe each quasar at roughly the 
same rest wavelength range, from Lya through C IV emis- 
sion. One exception is 0946+3009, which has a redshift too 
low for the Siiv emission to appear in our spectra. 

We supple ment our data with spectra from the SDSS 
Data Release 6 (|Adelman-McCarthv et alj|2008h for 8 of the 
quasars in our sample, for which the resolution is R ~ 2000 
(150 km s _1 ). These spectra cover the observed wavelength 
range 3800 to 9200 A, and we only include spectra that cover 
at least the Si IV through C IV emission. 

Table 1 summarizes the full dataset presented in this 
work, including the emission redshift, z em |j and the 'bal- 
nicity index' (BI) for each object (as calculated in Paper 
1). Any uncertainty in the redshift will not affect our com- 
parisons between Si IV and Civ or any of our other results. 



2 The values of z em were obtained from the NASA/IPAC Extra- 
galactic Database (NED), which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under con- 
tract with the National Aeronautics and Space Administration. 



The balnicity index, defined bv lWevmann et al.l |l99lj), is a 
measure of the strength of the BAL absorption and is calcu- 
lated as an EW in units of velocity. It quantifies blue-shifted 
C IV absorption between —25 000 and —3 000 km s~ that 
reaches at least 10 per cent below the continuum across a 
region at least 2000 km s _1 in width. The next four columns 
list the number of observations taken for each object at each 
observatory and then the total overall number of observa- 
tions. The final column lists the range in At covered for each 
quasar. 

Two of the quasars in our sample have BI=0, so they are 
not BAL quasars based on the balnicity index. They both 
contain broad absorption, but this absorption falls outside 
the velocity range, —25 000 to —3 000 km s~ , used to define 
BI. As noted in Paper 1 and discussed further in Section 3 
below, including these two objects in our sample does not 
affect any of our main results. 

2.2 Measuring BALs and their variability 

In Fig.[T] we plot spectra for all 24 objects, showing the long- 
term comparisons (At = 3.8—7.7 yr) between a Lick spec- 
trum and an MDM spectrum. The one exception is 2225- 
0534, for which we only have short-term Lick data. For each 
object, we plot the Civ absorption region in the top panel 
and the corresponding Si IV absorption region in the bottom 
panel. The velocity scale is based on the wavelengths of C IV 
and Si IV in the observed frame calculated from the redshifts 
given in Table 1. In order to compare the Civ and Siiv ab- 
sorption regions, we use the bluer line in both the C IV and 
Siiv doublet for the zero-points of the velocity scales, i.e., 
1548.20 A for Civ and 1393.76 A for Siiv. 

We adopt the velocity ranges over which C IV BAL ab- 
sorption occurs defined in Paper 1. These regions were de- 
fined based on the definition of BI, i.e. they must contain 
contiguous absorption that reaches JSlO per cent below the 
continuum across ^2000 km s — 1 . We apply the same defini- 
tion when defining the velocity ranges of Si IV BAL absorp- 
tion. 

In Paper 1, we defined a pseudo-continuum fit for the 
fiducial Lick observation used in the long-term analysis for 
each quasar by fitting a power- law to regions of the spectrum 
free of absorption and emission. The preferred spectral re- 
gions for the fits were 1270-1350 A and 1680-1800 A, but 
were adjusted to avoid emission and absorption features as 
much as possible or due to the limits of the wavelength cov- 
erage. We then fit the C IV emission lines, using between 1 
and 3 Gaussians to define the line profile. In Paper 1, we 
also fit the Si IV emission in cases where C IV absorption 
overlaps with the Si IV emission. For this work, we addition- 
ally fit the Si IV emission in the fiducial Lick observation for 
all the quasars. For the long-term comparisons below, we 
fit the Si IV emission for the MDM spectrum in cases where 
the emission line varied. Some special cases where there were 
difficulties with fitting the Si IV emission, such as 1011+0906 
and 1309—0536, are discussed further in Section 3.3 below. 

When comparing multiple epochs, we scaled all the 
spectra to the fiducial Lick spectrum used for the pseudo- 
continuum fit. We only fit the power-law continuum to one 
spectrum for each object, so any errors in this continuum 
fit will not effect our main variability results. To match the 
individual epochs for each quasar, we adopt a simple ver- 
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Figure 1. Spectra of all 24 quasars in our sample, showing the long-term comparisons (At = 3.8—7.7 yr) between a Lick Observatory 
spectrum (bold curves) and a recent MDM spectrum (thin curves). For each quasar, the Civ region is displayed in the top panel and the 
corresponding Si IV region is shown in the bottom panel. For 2225-0534, we only have short-term Lick data (see Table 1). The vertical 
flux scale applies to the Lick data, and the MDM spectrum has been scaled to match the Lick data in the continuum. The dashed curves 
show our pseudo-continuum fits. The horizontal bars indicate intervals of BAL absorption included in this study, and the shaded regions 
indicate intervals of variation within the BALs. We used binomial smoothing to improve the presentation of the spectra. The formal 
lcr errors are shown near the bottom of each panel. (The two variability intervals defined for 1524+5147 were labeled as one interval in 
Paper 1.) 
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Figure 1. continued. 
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Figure 1. continued. 
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Figure 1. continued. 
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Table 1. Quasar Data 



Name 


Z-em 


BI 


Lick 


SDSS 


MDM 


Total 


At 








1988-92 


2000-06 


2007-09 




(yrs) 


0019+0107 


2.130 


2290 


6 





1 


7 


0.08-5.79 


0043+0048 


2.137 


4330 


2 


2 


1 


5 


0.35-6.13 


0119+0310 


2.090 


6070 


2 





1 


3 


0.65-5.57 


0146+0142 


2.909 


5780 


2 





2 


1 


0.52-5.15 


0226-1024 


2.256 


7770 


1 





1 


2 


4.66 


0302+1705 


2.890 





2 





1 


3 


0.27-4.42 


0842+3431 


2.150 


4430 


6 


1 


3 


10 


0.06-5.87 


0846+1540 


2.928 





5 





3 


8 


0.04-4.93 


0903+1734 


2.771 


10700 


2 


1 


1 


7 


0.04-5.29 


0932+5006 


1.926 


7920 


•1 


1 


■1 


9 


0.05-6.98 


0946+3009 


1.221 


5550 


2 





3 


5 


0.11-8.16 


0957-0535 


1.810 


2670 


2 





2 


1 


0.11-6.21 


1011+0906 


2.268 


6100 


1 





3 


7 


0.10-5.94 


1232+1325 


2.364 


11000 


1 





2 


3 


0.35-5.93 


1246-0542 


2.236 


4810 


2 





2 


4 


0.40-5.90 


1303+3048 


1.770 


1390 


1 





1 


5 


0.05-6.10 


1309-0536 


2.224 


4690 


2 





2 


1 


0.68-6.19 


1331-0108 


1.876 


10400 


2 


1 


2 


5 


0.42-5.97 


1336+1335 


2.445 


7120 


1 





3 


4 


0.07-5.79 


1413+1143 


2.558 


6810 


2 


1 


2 


5 


0.26-5.61 


1423+5000 


2.252 


3060 


2 


1 


2 


5 


0.39-5.87 


1435+5005 


1.587 


11500 


2 





2 


4 


0.34-7.72 


1524+5147 


2.883 


1810 


3 


1 


3 


7 


0.04-5.14 


2225-0534 


1.981 


7920 


3 








3 


0.27-0.73 



tical scaling that matches the spectra along the continuum 
redwards of the C IV emission line (i.e. from 1560 A to the 
limit of the wavelength coverage) , between the Si IV and C IV 
emission (~1425-1515 A), and between the Lya and the Si IV 
emission (^1305-1315 A). For the few cases where a simple 
scaling did not produce a good match and there were dis- 
parities in the overall spectral shape between the compari- 
son spectra, we fit either a linear function (for 0903+1734, 
1413+1143, 1423+5000, and 1524+5147) or quadratic func- 
tion (for 0019+0107 and 1309-0536) to the ratio of the two 
spectra across regions that avoid the BALs. We then multi- 
plied this function by the SDSS or MDM spectrum to match 
the fiducial Lick data. 

With the spectra for each object matched, we used vi- 
sual inspection to identify velocity intervals with a width 
of at least 1200 km s _1 that varied. We identify intervals 
of variability separately for Civ and Siiv absorption. We 
then calculate the average flux and associated error for each 
candidate variable interval in each of the two epochs being 
compared. The error on the average flux is given by 



2 



t V"^ 2 

,J Z — / 



where ai is the error on an individual pixel, taken from the 
error arrays as displayed in Fig. 1, and n is the number of 
pixels in the candidate variable interval. We then calculate 
the flux difference between the two spectra and place an 
error on this flux difference using the error on the average 
flux from each epoch. We include all intervals of variability 
where the flux differences are at least 4a. Any interval which 
varied by at least 4a was readily identified by our initial 
inspection procedure. 

However, photon statistics alone are not sufficient for 



defining real variability, so we took a conservative approach, 
described in more detail in Paper 1, whereby we omit 
ambiguous cases of variability, even if they meet the 4a 
threshold. Flux calibrations, a poorly constrained contin- 
uum placement, and underlying emission- line variability can 
all add additional uncertainty to identifying variability. For 
example, in 1011+0906 and 1309-0536, there might be BAL 
absorption, and variability, on top of the Siiv emission line 
(see Fig. 1). However, it is too ambiguous to be included in 
this study. See Paper 1 for further examples of intervals of 
potential variability that were not included because of ad- 
ditional uncertainties and Section 3.3 below, where we com- 
ment further on certain individual quasars. We include nar- 
row intervals of variability such as the shaded region in Si IV 
in 0043+0048 and the shaded region in C IV in 1246-0542 
in Fig. 1, where the flux differences are 6.3er and 5.6<r, re- 
spectively. These regions meet the aforementioned thresh- 
olds, and the comparison spectra match well in regions of 
the continuum free of emission and absorption and on either 
side of the variability interval. We also include regions such 
as those shaded in 1011+0906 in Fig. 1 because even though 
the errors are slightly higher in the MDM spectra shown, the 
flux differences are still 7.7<r for the variable region in Civ 
and 6.3a and 7.9cr for the two regions of variability in Si IV. 
Overall our approach is designed to be conservative; we try 
to exclude marginal cases of variability to avoid overestimat- 
ing the true variability fractions. 

We calculated the absorption strength, A, of the BALs 
in our sample, where A is the fraction of the normalized 
continuum flux removed by absorption (0 ^ A ^ 1) within a 
specified velocity interval. These calculations are described 
in Paper 1. Very briefly, we divide each interval of variability 
and absorption, as defined above, into equal-sized bins of 
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Figure 2. The average normalized absorption strength, (A), in 
C IV versus the (A) in Si IV for each absorption bin in each quasar 
in the long-term subsample. The error bars represent the Ict errors 
based on photon statistics. 



width 1000 to 2000 km s~\ with the final bin size depending 
on the total velocity width of the specified interval. Then, for 
each quasar, we adopt the same bin size for the epochs being 
compared and calculate (^4) and AA in every individual bin. 

One of the difficulties in directly comparing C IV to Si IV 
absorption is the wider doublet separation in SilV (500 km 
s _1 in Civ versus 1900 km s _1 in Siiv). This can cause the 
Si IV absorption intervals to be wider than those in C IV. The 
only effect this should have on the variability results in Sec- 
tion 3 below is that there may be portions of Si IV absorption 
that are detected as variable but the corresponding velocity 
intervals in Civ may be too narrow to pass our variabil- 
ity threshold. We comment further on this in Section 3.1. 
We mark the regions defined as BAL absorption (horizon- 
tal bars) and variability (shaded rectangles) in Fig. [l] The 
absorption and variability regions for Civ were defined in 
Paper 1. We defined the absorption and variability regions 
for Si IV independently from what we found for C IV. See 
Section 3.1 for a full discussion of these results. 

In Fig. [2] we plot the relationship between absorption 
strength, (A), in Civ and (A) in SilV for the long-term 
spectra shown in Fig. 1. Each point represents a different 
absorption bin in an individual quasar, such that individual 
quasars contribute multiple points to this plot. The diago- 
nal line through the plot represents equal strength in both 
lines. Since we fix the velocity scale based on the bluer dou- 
blet member in C IV and in Si IV and the doublet separation 
is wider in Siiv, the edges of the Siiv absorption troughs 
tend to extend redward of the edges of the corresponding 
Civ troughs (see, for example, 0302+1705 in Fig. 1). As 
mentioned above, to calculate the (A) values, we divide each 
BAL into bins. We thus remove the redmost bin for each ab- 
sorption trough because Siiv may have a greater strength 
than C IV in those bins due to the wider doublet separation 
and not necessarily because the Siiv absorption is actually 



stronger than C IV. We also plot la error bars, calculated 
by using the error spectra shown in Fig. 1 and averaging 
over the velocity interval for each bin. We point out that 
non-statistical errors, e.g., from the continuum fitting, can 
increase the error in (A) measurements by up to 0.05 to 0.1. 
It is clear from Fig. [2] that absorption in C IV is roughly as 
strong or stronger than the corresponding Siiv absorption. 
In some cases, the C IV has no detectable corresponding Si IV 
absorption at all. 

There are almost no well-measured cases of C IV ab- 
sorption weaker than Siiv. We did find a few absorption 
bins where Siiv appeared to be stronger than Civ, but 
some of these bins might not actually have absorption due 
to Si IV, or due solely to Si IV. In about 10% of BAL spec- 
tra, there are lower ioniz ation lines, s uch as Cn A1335 and 
Aim A1855,1863 (|Trump et all 120061 ). The wavelength of 
Cn places any Cn absorption line at a velocity of < — 13 
500 km s _1 in Siiv velocity space. In order to look for in- 
terloping Cn lines, we checked all of our objects for another 
low ionization species, Al III, which exists redward of C IV, in 
a region of the spectrum relatively uncontaminated by other 
lines. We found that in two cases (1232+1325, 1331-0108), 
the velocity of the Al ill line puts the corresponding C II line 
within a Siiv BAL trough. We therefore removed the bins 
affected by C II from Fig. [21 and we exclude the contami- 
nated Siiv regions in these two quasars from the analysis 
below. There are still a few points below the one-to-one line 
in Fig. [2] However, these few points are mostly within 3<r of 
the one-to-one line and therefore are consistent with equal 
strength changes in C IV and Si IV. The one point that is just 
beyond 3<r from the line corresponds to the interval —3100 
to -1200 km s" 1 in 1303+3048. 



3 RESULTS 

3.1 Variability in Siiv versus Civ BALs 

In this section, we directly compare the variability of Siiv to 
Civ in the "long-term" dataset from Paper 1. This involves 
2 epochs of data for 23 quasars separated by 3.8 to 7.7 yrs. 
Our main goal is to discriminate between the possible causes 
of BAL variability. We begin by looking at what fractions of 
quasars exhibited C IV and Si IV variability to determine if 
Si IV varies more or less often than C IV. In Paper 1, we found 
a correlation between the incidence of C IV BAL variability 
and outflow velocity and absorption strength. Here we in- 
vestigate whether similar trends exists for Siiv BALs. We 
then look at the relationship between the incidence of C IV 
variability and Siiv strength. Last, we compare the change 
in strength for the two lines when they both vary. 

In the long-term dataset in Paper 1, 15 out of 23 quasars 
(65 per cent) exhibited Civ BAL variability and 11 out of 
19 (58 per cent) exhibited Si IV BAL variability, at any mea- 
sured velocity. We do not have data covering the Si IV region 
for 2 of our quasars, and another 2 quasars do not have Si IV 
BALs. This comparison between Civ and Siiv is compli- 
cated because the absorption in C IV is not always accompa- 
nied by corresponding absorption (e.g., at the same veloc- 
ities) in Siiv (Figs. 1 and 2). In addition, we are observa- 
tionally less sensitive to absorption and variability at high 
velocities in Si IV, compared to C IV, because those wave- 
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lengths can have poorer signal-to-noise ratios and larger un- 
certainties in the continuum placement caused by blends 
with underlying broad emission lines. Altogether this means 
we are more sensitive to variability in C IV than Si IV in our 
data set. 

To make a fair comparison between the incidence of 
variability in the two lines, we recalculate the above fractions 
while considering just the flow speeds at >— 20 000 km s _1 . 
We adopt this velocity as the cutoff because, in some spec- 
tra, there is emission due to O I at ^—20 500 km s" 1 in the 
Si IV absorption region (see also I Gibson et al.l 120101 ). With 
this additional restriction, we find that Si IV is more likely to 
vary than C IV. In particular, 35 per cent (8/23) of quasars 
exhibited Civ variability and 47 per cent (9/19) exhibited 
Si IV variability. The dramatic reduction in the C IV variabil- 
ity recorded this way, compared to the 65 per cent quoted 
above, is due to i) consideration of a narrower velocity range, 
and ii) the specific exclusion of high velocities, which are the 
most likely to show variability (Paper 1). Nearly half of the 
occurrences of C IV variability detected in our data set are at 
high velocities, i.e., v <— 20 000 km/s. The incidence of C IV 
variability further reduces to 31 per cent (6/19) if we only 
include the 19 quasars which have complete spectral cover- 
age across Siiv and have a Siiv BAL. The further decline 
in the C IV variability in this case probably occurs because 
the two quasars excluded for having no Si IV BAL have weak 
C IV lines, and weak C IV lines are more likely to vary than 
strong ones (Paper 1). Thus, we again removed Civ BALs 
that are more likely to vary. 

Overall, it is important to realize that a number of fac- 
tors can affect the measured incidence of BAL variability. 
Our comparisons show that, over matching velocity ranges, 
Si IV BALs have a significantly higher incidence of variabil- 
ity than C IV BALs. This difference is probably related to 
the different line strengths. In particular, the Siiv BALs are 
generally weaker than Civ BALs (Fig. 0, and weaker lines 
tend to be more variable (Paper 1 and Figs. [4] and [5] below). 

To more directly compare Siiv and Civ BAL variabil- 
ity, we examine the individual velocity intervals over which 
the variability occurs. We consider intervals at all velocities 
here and in the remainder of this section. Variability in C IV 
occurred in a total of 20 velocity intervals. Ten of these in- 
tervals have measurable Siiv absorption. Nine of these 10 
intervals, or 90 per cent, showed Siiv variations in the same 
sense (either getting stronger or weaker) as the C IV changes. 
There is only one interval (in 0119+0310) that exhibits vari- 
ability in C IV, without corresponding variability at the same 
velocities in Siiv. Conversely, we find long-term Siiv BAL 
variability in a total of 22 velocity intervals. All of the vari- 
able Siiv intervals have significant corresponding Civ ab- 
sorption, and 10 of these intervals showed C IV variations in 
the same sense as the Siiv (45 per cent). 

As mentioned in Section 2.2, Siiv has a wider doublet 
separation than Civ, causing some of the Siiv absorption 
and variability intervals to be wider than the corresponding 
C IV intervals. For the intervals of Si IV variability without 
corresponding C IV variability, we looked for any evidence of 
variability in C IV that was not included because the width 
of the candidate varying region was too narrow to meet our 
variability threshold (see Section 2.2). There is only one case 
where we detect a marginal narrow variability region in C IV 
corresponding to a region in Siiv classified as variable (in 
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Figure 3. The top two panels show the number of occurrences 
of Si IV (solid lines) and C IV (dashed lines) BAL absorption and 
variable absorption versus velocity. The third panel is the second 
panel divided by the first. The bottom panel shows the same curve 
from the third panel for Si IV with la error bars. 



0903+1734). Even if we were to count this as a variable C IV 
interval, still only 50 per cent of Si IV variability intervals 
would have corresponding C IV variability. Therefore, while 
91 per cent of the intervals of Civ variability had corre- 
sponding Si IV variability, only ~50 per cent of the intervals 
of Si IV variability had corresponding C IV variability. These 
results reinforce our nrain conclusion above, that Siiv BALs 
are more likely to vary than their C IV counterparts. 

Next we examine the dependence of Si IV variability on 
velocity and absorption strength, matching our analysis of 
Civ BALs in Paper 1. Figure 3 shows the incidence of Siiv 
absorption and Si IV variability versus velocity. For compar- 
ison, this figure also shows the corresponding data for C IV 
taken directly from fig. 3 in Paper 1 (dashed curves). The 
top two panels display the number of quasars with Si IV BAL 
absorption and with Siiv BAL variability at each velocity 
(solid curves). The third panel is the second panel divided by 
the top one, which gives the fraction of Si IV BALs that var- 
ied at each velocity. The top panel shows clearly that there 
are more C IV BALs than Si IV BALs at higher velocities. 

In Paper 1, we showed that the incidence of Civ vari- 
ability increases significantly with increasing velocity. This 
trend is not evide nt in the Si iv data. Figure 3 displays la er - 
ror bars based on lWilsonl l|l927fl and lAgresti fc Coulll l|l99ctl . 
These errors are based on counting statistics for the number 
of quasars with absorption and variability at each velocity. 
We perfornred a least-squares fit, and the slope of the Si IV 
data is —6.58 ± 4.33 x 10 -6 , in the formal unit, fraction per 
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Figure 4. The top two panels show the number of occurrences 
of Si IV BAL absorption and variability versus the average nor- 
malized absorption strength, (A), in Si IV. The third panel is the 
second panel divided by the first. The bottom panel shows the 
same curve from the third panel with ltr error bars overplotted. 
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The slope is non-zero at just a 1.5a significance. 



Therefore, while there might be a weak tendency for more 
variability in SilV at higher velocities (Figure 3), the trend 
is not statistically significant. 

To further match our analysis from Paper 1, we looked 
at the relationship between the incidence of Si IV variability 
and the absorption strength, (A), in SilV. As described in 
Paper 1 (and Section 2 above), we divide each BAL into 
~1200 km s _1 bins, and we treat each bin as an individ- 
ual occurrence of absorption. In Fig. 3] the top two panels 
show the number of these occurrences of Si IV absorption and 
the number of these occurrences that varied at each value 
of Si IV absorption strength, (A) . An individual quasar can 
contribute more than once to each point in the histogram. 
The third panel is the second panel divided by the top one, 
and the bottom panel shows the same curve from the third 
panel with error bars plotted, calculated in the same way as 
for Fig. [3] We find only a weak trend between Si IV variabil- 
ity and Si IV strength. The slope of the plot is -0.315±0.090 
fraction per unit absorption strength, which is non-zero at a 
3.5<t significance. This is much weaker than the trend be- 
tween the incidence of C IV variability and C IV strength 
found in Paper 1. We overplot this curve for Civ from fig. 
5 of Paper 1 in the third panel of Fig. 4 (dashed curve). 
This indicates that the occurrence of variability in SilV is 
less sensitive to the strength of the line than the occurrence 
of variability in C IV. 

Next, we looked at the relationship between the inci- 
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Figure 5. The top two panels show the number of occurrences 
of C IV BAL absorption and variability versus the average nor- 
malized absorption strength, (A), at the same velocities in Si IV. 
The third panel is the second panel divided by the first, and the 
bottom panel shows the same curve from the third panel with la 
error bars. In the third panel, we overplot the fraction of occur- 
rences of C IV absorption that varied versus (A) in C IV. 



dence of C IV BAL variability and the strength of the Si IV 
absorption, (A), at the same velocities. Si is known to be 
less abundant than C, when solar abundances are assumed, 
and the high ionization typical in BALs favors C IV (see Sec- 
tion 5). Therefore, the optical depth in Civ is higher than 
in Si IV, so the stronger the SilV absorption is, the more 
likely C IV is to be saturated. In Fig. [5] we plot the num- 
ber of occurrences of C IV absorption that occur at the same 
velocities in the same spectra as each SilV (A) value and 
then the number of these occurrences that varied in the sec- 
ond panel. As in Fig. 01 an individual quasar can contribute 
more than once to each point in the histogram. The third 
panel is the second panel divided by the top one, showing 
the fraction of occurrences of C IV absorption that varied at 
each Si IV absorption strength value. As in Figs. [3] and U the 
bottom panel of Fig. [5] shows the la error bars. The slope 
of these points is —0.643 ± 0.066 fraction per unit absorp- 
tion strength, which is non-zero at a 10a significance. Fig. [5] 
thus indicates that the incidence of C IV variability decreases 
with increasing SilV absorption strength. Therefore, when 
the Si IV absorption is stronger and the C IV absorption is 
more likely to be saturated, the incidence of C IV variability 
decreases. In fact, whenever the SilV absorption strength is 
greater than 0.5, the corresponding Civ absorption at the 
same velocity does not vary. 

As a further comparison to Paper 1, we overplot in the 
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Figure 6. The change in strength of Civ BAL absorption versus 
the change in strength of Si IV BAL absorption in velocity inter- 
vals where both lines varied. The error bars are calculated as in 
Fig.E 
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Figure 7. The fractional change in strength of C IV BAL absorp- 
tion versus the fractional change in strength of Si IV BAL absorp- 
tion in velocity intervals where both lines varied. The error bars 
are calculated as in Fig. [2] 



third panel of Fig. [5] the fraction of occurrences of C IV ab- 
sorption that varied versus C IV absorption strength (dashed 
curve). In Paper 1, we concluded that, for Civ BALs, weaker 
lines are more likely to vary than stronger lines. Fig.[5]shows 
that C IV lines are even more likely to vary when the corre- 
sponding Si IV lines are also weak. 

We also investigate the change in strength, | Aj4|, in C IV 
compared to the corresponding changes in SilV in the same 
velocity interval (Fig. [6]). There are 6 quasars in which the 
velocity intervals of C IV and SilV variability overlap, and we 
only compare the velocity intervals where both lines varied. 
Each point represents one bin in one of these quasars, as de- 
scribed above for Fig. [5] The diagonal line through this plot 
corresponds to equal strength changes in both lines. There is 
clearly no correlation between the strength changes in C IV 
versus SilV evident in this figure. Despite SilV varying more 
often than C IV, the strength changes in Si IV are not always 
greater than in C IV. 

Finally, Fig. [7] shows the fractional change in strength, 
|AA|/(yl), in Civ compared to |AA|/(j4) in Siiv, again for 
velocity intervals where both lines varied. As in Fig. [U there 
are points lying above the line, representing greater strength 
changes in C IV than in Si IV. However, there is a weak trend 
towards greater fractional change in strength in Si IV, which 
is consistent with the correlation found between fract ional 
change in EW in Civ and SilV in lGibson et al.l lJ2010h . 

As mentioned in Section 2.1, we include two quasars 
that have BI=0 because they do have broad absorption, but 
it falls outside the velocity range, —25 000 to —3 000 km s , 
used in the strict definition of BI. We also include broad ab- 
sorption in other quasars in our sample that falls outside 
of this velocity range. Their inclusion has minimal impact 
on our results because most of the Si IV broad absorption in 
our sample falls within this velocity range. For the quasars 
with BI=0, 0846+1540 does not contain any Siiv broad ab- 



sorption at all, and 0302+1705 contains broad absorption at 
low- velocity, which did not vary in either C IV or Si iv. 

We can summarize our comparisons between the C IV 
and Siiv BAL variabilities as follows: 1) Siiv BALs are more 
likely to vary than C IV BALs. The fractions of quasars show- 
ing variability in our long-term 2-epoch data set are 31% in 
C IV and 47% in Si IV if we consider only the well-measured 
velocity range v > —20 000 km/s and include only quasars 
with both Si IV and C IV BALs detected. 2) The variabilities 
usually occur in just portions of the BAL troughs. 3) When 
changes occur in both Si IV and C IV, they always occur in 
the same sense (i.e., with both lines getting either weaker or 
stronger). They also occur in overlapping but not necessarily 
identical velocity ranges. 4) The trend for a higher incidence 
of C IV variability at higher velocities, which we reported in 
Paper 1, is not clearly evident in the Siiv data. Finally, 5) 
there is no correlation between absorption strength changes 
in Civ versus Siiv when they both vary; although, there is 
a weak trend towards greater fractional change in strength 
in Siiv. 



3.2 Multi-Epoch Monitoring of BALQSOs 

We now expand our analysis to the full dataset, which in- 
cludes 2 to 10 epochs of data for each object (Table 1). 
Including all of these epochs and considering all measured 
velocities, the fraction of quasars that showed C IV BAL vari- 
ability is 83 per cent. This is a significant increase from 
the 65 per cent we derived considering only two long-term 
epochs, or the 39 per cent derived from only two short- 
term epochs (Section 3.1 and Paper 1). Clearly, including 
more epochs of data increases the observed variability frac- 
tions. Moreover, these larger variability fractions apply to 
roughly the same time-frame as our 2-epoch long-term data 
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set. Therefore, the multi-epoch data did not find new oc- 
currences of variability at some other time; they identified 
variability missed by the 2-epoch measurements. 

To investigate the multi-epoch behaviours of BAL vari- 
ability, we compared all the spectra obtained for each 
quasar. From one object to another, there are large differ- 
ences in the widths of the varying regions and the amplitudes 
of the changes (e.g., see Figs. 1 and [6]). However, there are 
certain general trends that most, if not all, the quasars fol- 
low. In particular, the variability almost always occurred 
within just a portion of a BAL and not in the entire trough. 
In nearly all the quasars, the variability occurred over the 
same velocity interval(s) between each epoch. Finally, when 
there are multiple velocity intervals of variability within the 
same quasar, the changes in these separate intervals almost 
always occur in the same sense. Similarly, as described in 
Section 3.1, when there is variability in both Civ and Siiv, 
they also vary in the same sense. 

We also find no clear evidence for velocity shifts that 
would be indicative of acceleration or deceleration in the 
flows. The constraints on velocity shifts are difficult to quan- 
tify in BALs because there can be complex profile variabil- 
ities, but we specifically search for and did not find cases 
where a distinct absorption feature preserved its identity 
while shifting in velocity. Despite the large outflow veloc- 
ities, there is no clear evidence to date for acceleration or 
deceleration in BALs , or in any other outflow lines (i .e. NALs 
and mini-BALs; e.g. iRodrfguez Hidalgo et al.ll201ll ). 

We highlight below a few well-sampled cases to illus- 
trate these general trends in the data. Fig. [S] shows specif- 
ically the quasars for which we have at least 7 measured 
epochs including one from the SDSS, which helps to span the 
time gap between the early Lick data and our recent MDM 
observations. For each object, the top panel shows the Civ 
BAL(s) and the bottom panel shows the Siiv BAL(s). The 
blue curves show the early Lick data, the red curves show 
the intermediate SDSS data, and the green curves show the 
MDM spectra. We note that in 1524+5147 there is an O I 
emission line centered at ~— 20 500 km s _1 in the Si IV panel. 

The bold bars in Fig. [5] mark intervals of variability 
identified for C IV. The velocity ranges are guided by the 
intervals defined for C IV in Paper 1 and are adjusted to 
cover the core of the varying region and avoid the edges 
where the variability is less pronounced. These same veloc- 
ity ranges are applied to all the epochs plotted here and 
to Siiv for comparison. The one exception is the thin bar 
marking a region of variability in Siiv in 0842+3431 with 
minimal corresponding variability in C IV. For all of these 
varying regions, we calculate the absorption strength, A, for 
the defined velocity intervals in each epoch and then plot 
these A values versus time in Fig. [9] 

We plot A instead of EW, in Fig. [9l in order to high- 
light the intervals that varied. Using EW would dilute these 
changes in strength. The different colors correspond to dif- 
ferent velocity intervals. The dashed and dotted lines rep- 
resent changes in Civ and Siiv, respectively. We note that 
these lines do not represent how the A values changed be- 
tween epochs. They simply connect the measurements from 
different epochs in order to aide the eye. 

In Fig. [51 the spectra of 0842+3431 (in the Civ BAL), 
0903+1734, and 1524+5147 show clearly how portions of 
BALs can vary. In 0842+3431, there is significant variabil- 



ity in both the blue side and the red side of the Civ BAL; 
although in Si IV, the entire BAL varies (see Section ^. 3. 3l be- 
low). In 0903+1734, there is significant variability at higher 
outflow velocities, but at lower velocities, there is no vari- 
ability. This is consistent with the result from Paper 1 that 
there is a higher incidence of variability at higher velocities. 
The variable regions in 1524+5147 cover most, but not all, 
of the BAL. In contrast to the general trends in Paper 1, the 
blue-most portion of the BAL did not vary. In 0932+5006, 
however, the entire C IV BALs vary at the highest velocities. 
In terms of Civ to Siiv comparisons, 1524+5147 is a case 
where there is weak corresponding absorption, and variabil- 
ity, in Si iv, but no Si iv BAL. And, 0932+5006 shows clearly 
a case where a Siiv BAL varied, but the Civ BAL did not. 
The plots in Fig. U for 0842+3431, 0903+1734, and 
0932+5006 show how the change in A is not always mono- 
tonic. The same BAL can grow deeper from one epoch to an- 
other, then b ecome shallower aga in. This is consistent with 
the results of iGibson et al.l (J2010r) . Furthermore, the change 
in the A value from one epoch to another generally occurs 
in the same direction (either positive or negative) in both 
C iv and Si iv, which is consistent with the results of Section 
3.1. In 0842+3431, 0903+1734, and 1524+5147, where there 
are two separate intervals of C iv variability, the change in 
strength occurs in the same direction for both intervals. The 
high-velocity BALs in 0932+5006 vary in concert starting 
with the 1989.84 epoch through 2008.03. At the earliest and 
latest epochs, they do not clearly vary in concert, but the 
changes in A are small and could be affected by changes in 
the underlying Siiv emission line. We comment further on 
these BALs in 0932+5006 in Section |3~3H 



3.3 Notes on Individual Quasars 

In this section, we comment on individual quasars that are 
cases of special scientific interest. We also comment on cases 
where there were specific issues in the analysis or measure- 
ments that result in larger uncertainties. 



3.3.1 0119+0310 

0119+0310 is the only quasar for which we record Civ 
BAL variations without corresponding changes in Siiv in 
our long-term sample (Fig. 1 and Section 3.1). However, 
these results are very tentative because the Si iv absorption 
is poorly measured across the velocities that varied in C IV. 
The two long-term spectra for this object, plotted in Fig. 1, 
have a lower signal-to-noise level than most of the other data 
in our sample. Furthermore, the Siiv absorption at the ve- 
locity of Civ variability (~— 7500 km s _1 ) is very weak, and 
if the continuum fit is off by even ~5 per cent, this region in 
Siiv might not be considered part of the Siiv BAL. If this 
region is not part of the BAL, then we would not include 
it in the comparison of C iv to Si iv. Therefore, while we 
have several well-measured cases of Si iv variability without 
corresponding C iv variability, we only have this one poorly- 
measured case of C IV variability with no corresponding Si iv 
variability. 

This quasar also appears to differ from most of the other 
quasars in the sample in that the different regions of C iv 
variability do not vary in the same sense (Fig. 1). The two 
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Figure 8. Spectra of the C IV (top panel) and Si IV (bottom panel) BALs in four well-sampled quasars from our sample, after smoothing 
three times with a binomial function. The blue curves are the Lick spectra, red curves are SDSS spectra, and green curves are MDM 
spectra. The bold bars mark varying regions identified for C IV. The average error for each spectrum is shown in the top panel for each 
quasar, where the height of the error bar represents ±1<T. 
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Figure 9. The absorption strength in both Civ (dashed lines) and Si IV (dotted lines) as a function of time in the velocity intervals 
indicated by bars in Fig. [8] 



higher velocity variable regions both increase in strength 
between the Lick and MDM observations, while the lower 
velocity variable region decreases in strength. As mentioned 
above, this is one of our least well-measured quasars, so this 
is a tentative result. We find just two other cases where 
two regions of C IV variability vary in opposite directions 
(0146+0142 and 1423+5000). 



3.3.2 0146+0142 

We first note that this object has a high-velocity Civ BAL, 
so the BAL that appears just redward of the marked Si IV 
BAL in Fig. 1 is actually Civ. As mentioned in Paper 1, 
we can confirm that this is high- velocity C IV, and not Si IV, 
because if it were Si IV absorption, we should see correspond- 



ing lo w- velocity C IV absorption (see also fig. l. lKorista et al.l 
1993). We also have further confirmation that this is high- 
velocity C IV absorption because we find evidence of corre- 
sponding high- velocity SilV absorption on top of the Lycv 
emission line. Fig. [TO] shows the two long-term spectra for 
0146+0142 with the two right-most shaded regions mark- 
ing the C IV variability and the two left-most shaded regions 
marking the corresponding velocities (but not necessarily 
the entire variable regions) in Si IV. 

Another interesting note about 0146+0142 is that in our 
short-term data, there are two separate regions of C IV BAL 
variability, but they do not vary in the same sense. In Fig. 1111 
we plot the two short-term epochs with these two variable 
regions shaded. The redder varying interval at ~ — 28 300 km 
s _1 increases in strength, while the bluer interval at ~— 36 
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Figure 10. The two long-term epochs for 0146+0142, showing 
the full spectrum from the Lycr emission through the C IV emis- 
sion. The shading here differs from Fig. 1, with the right-most 
shaded regions marking the C IV variability and the left-most 
shaded regions showing the corresponding velocities in Si IV (but 
not the exact regions of Si IV variability) . This figure shows evi- 
dence of Si IV BAL variability on top of the Lya emission line. 
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Figure 11. The two short-term epochs for 0146+0142, with the 
two shaded regions marking C IV variability. The absorption in 
these two regions vary in opposite directions, with one region 
becoming weaker while the other becomes stronger. 



600 km s 1 decreases in strength. As mentioned in Section 
13.3.11 this is one of just three cases showing this behavior. 



3.3.3 0842+3431 

In 0842+3431, there is significant variability in two distinct 
regions of the C IV BAL trough (marked by bold bars in 
Fig. |8}, while the entire Siiv trough varies. However, the 
bluer, and more variable, region in the C IV trough (the left- 
most bold bar in Fig. [H| corresponds to a region of weak 
absorption and variability in Si IV. Some of the variability in 
the redmost portion of the Si IV trough (the right-most bold 
bar in Fig. [5J may be connected with the variability in the 
red side of the C IV BAL trough, but the variability in the 
core of the Si IV trough, marked by the thin bar in Fig. [HI 



cannot be explained by the wider Si IV doublet separation 
alone. Some of the Siiv variability is occurring at different 
velocities than the C IV variability. Nonetheless, as seen in 
Fig. [9l the changes in strength in all three marked regions 
of the Siiv trough occur in the same sense as the changes 
in strength in the two distinct varying regions in C IV. 

Another interesting note about 0842+3431 is that we 
identified it as variable in the short-term, but not in the long- 
term, in Paper 1. For the long-term comparison in Paper 1, 
we used the 1990.90 and the 2008.35 observations. However, 
Figs. [8] and [9] clearly show that the C IV BAL varied. The 
strengh of the BAL is weaker in 2007.04 than in 1990.90, 
but the BAL becomes stronger again by 2008.35. Therefore, 
when looking at just the 1990.90 and 2008.35 observations, 
it appears as if the BAL did not vary at all. This shows 
how variable profiles can return to a previous state and that 
variable BALs can be missed in 2-epoch studies. 



3.3.4 0932+5006 

As in 0146+0142, we note that the absorption trough that 
appears just redward of the Si IV BAL in 0932+5006 in Fig. 
1 is actually a high- velocity Civ BAL (see Fig. [8j. 

In 0932+5006, there is C IV absorption overlapping 
Si IV emission. The variability in these two detectable BAL 
troughs has a different behaviour than the variability in the 
other quasars in our sample. When looking at the spectra 
(Fig. [8]), the two high- velocity BALs in the 2003.01 spectrum 
appear to be offset in velocity from the BALs in the other 
epochs. While this could be indicative of a shift in velocity 
of the BALs, this apparent offset could also be due to part 
of the trough weakening while the other part strengthens. 
This complicates the measurement of A for Fig. [5] because 
any measurement of A in a fixed velocity interval for each of 
these two high-velocity BALs does not accurately represent 
how the line changed. Furthermore, the Siiv emission line 
itself could be variable, which complicates any analysis of 
C IV BAL variability in this velocity range. 



3.3.5 0946+3009 

0946+3009 is the one object in our sample with a redshift 
too low for our MDM spectra to cover the entire C IV region 
out to the Siiv emission line. The spectra only go as blue 
as ~ — 19000 km s _1 . The absorption and variability regions 
that we define for this quasar do not extend all the way to 
the edge of the MDM spectrum in order to avoid any un- 
certainties there. The detection of variability in this quasar 
is secure because we have an additional MDM spectrum of 
this quasar that matches the MDM spectrum shown in Fig. 
1. 



3.3.6 1011+0906 

As mentioned in Section 2.2, some of the quasars in our sam- 
ple have low-ionization BALs. We searched for Al ill lines in 
our sample and used the velocity of the Al ill line to deter- 
mine the location of Cn. 1011+0906 has an Aim BAL, but 
the velocity of the absorption puts C II blueward of the Si IV 
absorption. Therefore, if there is any C II absorption in this 
object, it does not affect our measurements of the Siiv BAL. 
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The Si IV broad emission line (BEL) in this quasar is 
mostly absorbed by the Civ BAL. We fit the Siiv BEL 
using the procedure defined in Paper 1 for cases like this. We 
take the Civ fit, increase the FWHM based on the greater 
doublet separation in Siiv, and place it at the wavelength 
where the Siiv emission should be. However, there is still 
some slight emission blueward of the Siiv BEL fit. This 
extra emission may be part of the wing of the Siiv BEL, 
or the underlying power-law continuum fit might be slightly 
too low. However, the Siiv BAL is located at a high enough 
velocity that any error in the Siiv BEL fit should have a 
negligible effect on our measurements of the BAL and its 
variability. 



3.3.10 1331-0108 

As in 1232+1325, 1331-0108 has an Aim BAL at a velocity 
that places the corresponding C II absorption within the Si IV 
BAL. We therefore omit the velocity range —23 800 to —15 
300 km s _1 in the Siiv region from our analysis. 

While analyzing the Siiv region in 1331—0108, we no- 
ticed that the pseudo-continuum fit defined in Paper 1 
needed to adjusted. Like 1309—0536, 1331—0108 has very 
broad BALs, which makes fitting a continuum difficult. The 
slope of the fit for 1331—0108 is now slightly steeper than 
the fit used in Paper 1, increasing the measured A values 
for C IV on average by just ~7 per cent throughout most of 
the trough and up to ~30 per cent at the highest velocities, 
where the absorption is much weaker. 



3.3.7 1232+1325 

1232+1325 has an Al ill BAL which puts C II within the Si IV 
region. The C II BAL is in the velocity range —25 500 to —19 
600 km s _1 (see also Fig. 1), and we omit these velocities 
from our analysis. 



3.3.8 1303+3048 

1303+3048 is a BAL quasar that also contained a C IV mini- 
BAL at ~ — 18 500 km s _1 when first observed at Lick (Fig. 
1). We only have one Lick observation of this object, but the 
MDM observations show that the mini-BAL widened and 
increased in strength to become a BAL. A BAL emerges at 
the same velocities in Siiv as well. The lower- velocity BAL 
in 1303+3048 is visible in C iv in the Lick data, but appears 
as only weak absorption in Si IV. However, between the Lick 
and MDM epochs, a Si iv BAL emerges and the variability in 
the Si iv absorption extends to lower velocities than the C IV 
variability. The variability at all velocities in this quasar in 
both C iv and Si iv occurs in the same sense; the absorption 
increases in strength. 



3.3.9 1309-0536 

As in 1011+0906, the Siiv BEL in 1309-0536 is heavily ab- 
sorbed by a C iv BAL. We used the same procedure that 
we used for 1011+0906 to fit the Siiv BEL, and we found 
what appeared to be significant emission blueward of the 
Siiv BEL fit. In this case, the underlying powerlaw contin- 
uum fit did not have the correct slope, so we made a slight 
adjustment to the continuum fit. Adjusting the powerlaw 
continuum fit caused on average an increase in A of ~6 per 
cent throughout most of the C iv trough, compared to the 
measurements in Paper 1. The new Si iv BEL fit increased A 
at the highest velocities in Civ by up to a factor of 2. Even 
with this adjustment, there still appears to be some slight 
emission blueward of the Siiv BEL fit, but, as in 1011+0906, 
the Siiv BAL in 1309—0536 is at a high enough velocity 
that errors in the emission fit should not have much effect 
on measurements of the BAL. This quasar also did not vary 
in either Si iv or C iv, so any measurement errors for this 
quasar will not effect any of our results comparing Si iv and 
Civ variability properties (e.g., Figs. [6] or [7|. 



3.3.11 1423+5000 

1423+5000 is another quasar where there were two C iv 
BALs that varied, but one BAL increased in strength, while 
the other weakened. This quasar varied between the Lick and 
SDSS observations, but we did not detect any variability in 
the long-term analysis in Paper 1. 0119+0310, 0146+0142, 
and 1423+5000 are the only quasars in our sample where we 
see two separate varying regions in C iv that did not vary in 
the same sense. 



3.3.12 1435+5005 

1435+5005 has Aim absorption that is either a weak BAL 
or strong mini-BAL. However, the velocity of the absorption 
places Cn blueward of the Siiv absorption. 

We also note that the signal-to-noise level in the data 
for 1435+5005 decreases rapidly at bluer wavelengths. We 
therefore we do not include the spectral region blueward of 
— 12 100 km s" 1 in Siiv in our analysis in Section 3.1. 



4 SUMMARY OF RESULTS 

This is the second paper in a 3-part series to analyze the 
BAL variabilities in a sample of 24 BAL quasars measured 
originally by Barlow (1993) at the Lick Observatory in 1988- 
1992. We supplement those data with spectra from the SDSS 
archives (for 8 quasars) and our own measurements obtained 
at the MDM observatory (Table 1). In Paper 1 we discussed 
the variability properties of C iv A1549 measured in just 
two epochs that span a "short-term" (0.35—0.75 yr) and 
a "long-term" (3.8—7.7 yrs) time interval. Here we build 
upon that work by including our full multi-epoch data set 
for these same quasars and making detailed comparisons 
between the Si iv and C iv BAL behaviors. Our main results 
are the following: 

(1) BAL variability usually occurred in only portions of 
the BAL troughs (Paper 1; Section 3.3.2). 

(2) In the long-term interval, 65 per cent of the BAL 
quasars in our sample showed C iv BAL variability while 
only 39 per cent varied in the short-term (Paper 1). 

(3) C iv variability occurs more often at higher velocities 
and in shallower absorption troughs (or shallower portions 
of absorption troughs; Paper 1). 
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(4) In rare cases, BAL features appear, disappear, or 
change to or from narrower mini-BAL features (Paper 1; 
Section EHTH)) . 

(5) C IV BALs in our data are as strong or stronger than 
Si IV BALs at all velocities (in all well-measured cases; Fig. 

(6) Si IV BALs are more likely to vary than Civ BALs. 
For example, when looking at flow speeds >— 20 000 km s , 
47 per cent of the quasars in our sample exhibited Si IV vari- 
ability while 31 per cent exhibited C IV variability (Section 
3.1). The greater variability in Siiv is likely due to a com- 
bination of items (3) and (5) above; weaker lines are more 
likely to vary, and Si iv tends to be weaker than C iv. 

(7) Variability in Siiv can occur without corresponding 
changes in C iv at the same velocities. ~50 per cent of the 
variable Si iv regions did not have corresponding C iv vari- 
ability at the same velocities. However, in only one poorly- 
measured case were changes in C iv not matched by Si IV 
(Sections 3.1 and 13.3. 1) , 

(8) At BAL velocities where both Civ and Siiv varied, 
the changes always occurred in the same sense (Section 3.1). 

(9) We do not find any correlation between the absolute 
change in strength in C iv and in Si iv (Fig. |S} , but the 
fractional change in strength tends to be greater in Siiv 
than in Civ (Fig. 0. 

(10) When additional observing epochs are included (e.g., 
our full data set; Section 3.2), the fraction of C iv BALs that 
varied at any velocity increases from 65 per cent to 83 per 
cent. This increase was caused by variations missed in the 
2-epoch comparisons in Paper 1. 

(11) BAL changes at different velocities in the same ion 
almost always occurred in the same sense (getting weaker 
or stronger) but not generally by the same amount (Section 
3.2). We find just 3 cases that show evidence for one Civ 
BAL weakening while another strengthens within the same 
object (Sections I337T1 l3~3"l?l and l3.3.1ip . 

(12) The multi-epoch data also show that the BAL 
changes across 0.04—8.2 years in the rest frame were not 
generally monotonic (Section 3.2). Thus, the characteristic 
time-scale for significant line variations, and (perhaps) for 
structural changes in the outflows, is less than a few years. 

(13) With more epochs added, we still do not find clear 
evidence for acceleration or deceleration in the BAL outflows 
(Section 3.2). 



5 DISCUSSION 

The BAL variability data provide important constraints on 
the outflow physical properties. However, the information 
we derive depends critically on what causes the BAL varia- 
tions. In this section we discuss pros and cons of two compet- 
ing scenarios, namely, 1) fluctuations in the far-UV contin- 
uum flux that cause global changes in the outflow ionization, 
and 2) outflow clouds moving across our lines-of-sight to the 
quasar continuum source. 

An important part of this discussion is the BAL optical 
depths, which can be much larger than they appear in the 
spectrum if t he absorbers cover only part of the b ackground 
light source IjHaman n 1998; Ham ann et al.l l2008). Compar- 
isons between the C iv A1549 and Si iv A1400 BALs can help 
because these lines probe slightly different ionizations with 



potentially very different line optical depths. For example, 
in a simple situation with solar abundances and an ion ra- 
tio equal to the abundance ratio, i.e., Siiv/Civ = Si/C, 
the optical depth in Siiv A1400 would be ~3.4 times less 
than C IV A 1549 dHamannl ll997J ; lHamann fc Ferlandl Il999l ; 
lAsplund et al.l|2009l ). In actual BAL flows, the relative Siiv 
optical depth should be even lower because BAL ionization 
tends to be high and thus favors Civ. We cannot make 
specific comparisons without specific knowledge of the ab- 
sorber ionizations. However, if we reasonably assume that 
the ionization is at least as high as that needed for a maxi- 
mum Civ/C ratio (e.g., in a gas that is photoionized by the 
quasar and optically thin in the Lyman continuum - fig. Al 
in Hamann et al. 2011), then the Siiv optical depths should 
be >8 times smaller than C iv. 



5.1 Changing Ionization 

When there is variability in different velocity intervals within 
the same BAL or within multiple BALs in the same quasar, 
the changes almost always occur in the same sense (e.g., 
0842+3431 and 0903+1734 - Figs. [S] and 0. Studies of nar- 
row absorption line (NAL) variability have o bserved multi- 
ple NALs in a given quasar varying in concert (JMisawa et al.l 
120071 ; lHamann et al.ll201ll '). lHamann et al.l J201ll ) found co- 
ordinated line variations in five NAL systems in a single 
quasar. They argue that the most likely explanation for 
this is a global change in ionization. If there are changes in 
the ionizing flux incident on the entire outflow, then global 
changes in ionization should occur. While the connection be- 
tween NALs and BALs is unclear, this argument can be ap- 
plied to BALs as well. Absorbing regions at different veloc- 
ities have different radial distances from the central SMBH. 
They are therefore spatially distinct, even if they are part 
of the same larger outflow structure. A change in covering 
fraction due to moving clouds is unlikely in cases such as 
0842+3431 and 0903+1734 because it would require coor- 
dinated movements among multiple absorbing structures at 
different outflow velocities and radii. 

To further investigate this scenario, for simplicity, we 
consider a system with homogeneous clouds outflowing from 
the accretion disk, with no transverse motion across our line- 
of-sight to the quasar. A change in ionization will cause the 
optical depths in the lines to change. As mentioned above, 
the optical depths in C iv are higher than in Si iv, so Si IV 
would be more susceptible to changes in ionization. There- 
fore, it is more likely for Si iv to vary than C iv in this sce- 
nario because C iv is more likely to be saturated. This gen- 
erally matches our results since we find that Siiv is more 
variable than C iv. We find only one case of C iv variability 
unaccompanied by Si iv variability, and it is a very tentative 
result (Section 13.3. l|l . 

This scenario becomes a little more complicated when 
considering that typically variability only occurs in portions 
of BAL troughs, r ather than entire BAL troughs varying 
II Gibson et al. 2008; Paper 1). As mentioned above, a change 
in ionization should cause more global changes, rather than 
changes in small, discrete velocity intervals. It is possible 
that the variable regions in the troughs have moderate or 
low optical depths, while the non-variable sections are too 
saturated to respond to modest changes in the ionization 
and line optical depths. We have evidence from Paper 1 
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and Fig. [5] that weaker lines, or weaker portions of lines, 
are more likely to vary. However, we also found in Paper 1 
that variability is more common at higher velocities, where 
the absorption tends to be weaker, so it is difficult to say 
whether it is the higher velocity or the weaker absorption 
strength that is the root cause of the variability. If it is 
true that weaker portions of lines, which are least likely to 
be saturated, are more likely to vary, regardless of outflow 
velocity, then this would support the changing ionization 
scenario. 

However, if it is true that weaker portions of lines are 
less saturated and thus more likely to vary, it is unclear 
why, for example, the weak blue wing of the BAL trough in 
0842+3431 does not vary. If changing ionization is causing 
the variability in this quasar, then the wings of the line must 
be saturated while the portions of the line adjacent to the 
wings are not saturated. There are other examples of similar 
behavior. In 1524+5147, the strongest variability occurs in 
the deepest segment of the BAL, and there is weak or no 
variability in the weakest segments of the trough (Figs. 1 
and[8]). Similarly, the variability in 1011+0906 occurred near 
the core of the line, while the wings did not vary (Fig. 1). 

In order for changing ionization to cause variability in 
just the deepest portions of BAL troughs, as in the exam- 
ples given above, there must be velocity-dependent cover- 
ing fractions with velocity-dependent optical depths. In this 
way, even the weak wings can be highly saturated. There 
is evidence in the literature that both optical depth and 
covering fractions can have complex vel ocity-de pendent be- 
haviors (Barlow & Sargent 1997; Hama nn et al.lll997l. 



Ganguly et al J 1 19991: Ide Kool et al][2002l ; iGabel et all 
Gabel et al.ll2006l : lArav et al.ll200; 



2001; 



2005, 



The true optical depths and covering fractions are diffi- 
cult to measure for BALs. In our data, 1413+1143 provides 
direct evidence for velocity-dependent optical depths if the 
line variations are caused by ionization changes. At the core 
of the Si IV trough there are two dips at the Siiv doublet 
separation (Figure 1). The doublet ratio is roughly one-to- 
one, indicating saturation at these velocities and therefore 
little or no sensitivity to changes in continuum flux. This 
part of the trough did not vary. However, this saturated 
doublet is surrounded by variability at higher and lower ve- 
locities in Si IV. In C IV, which should have generally larger 
optical depths, variability occurs only at higher velocities. 
This behavior is at least suggestive of lower optical depths 
(non-saturated absorption) at the variable velocities. 

One further piece of evidence for the changing ioniza- 
tion scenario comes from the multi-epoch data in Section 
3.2, which show that changes in BAL strength ar e not nec- 
essarily monotonic (see also, iGibson et al.l 120101 ). In fact, 
in 0842+3431, the absorption trough varied, and then in 
our last MDM observation it returned to the same strength 
it had in one of the first Lick observations. In order for a 
change in covering fraction to have caused the variability in 
0842+3431, the cloud movements would have to be repeat- 
able, in addition to being coordinated at different velocities 
corresponding to different spatial locations. A change in ion- 
ization is a more likely explanation because continu um flux 
varia tions are not necessarily monotonic either (e.g. iBarlowl 
1 19931 ). 

If a change in ionization does indeed cause the variabil- 
ity we detect, then there should be a connection between 



changes in continuum flux and BAL variability. However, the 
results from previous studies have been mixed. IBarlowl 1 19931 
found some evidence for a correlation between continuum 
variability and BAL variations, at least in certain individual 
quasa rs, while other studies have not found a strong correla- 
tion JBarlow et al.lll992l iLundgren et al.ll2007l : IGibson et all 
2008). However, all of these studies look at near-UV flux 
variations and little is known about the far-UV variability 
properties of luminous quasars. It is the far-UV flux that 
is the source of the ionizing radiation. Therefore, these re- 
sults do not rule out ionization changes as a cause of BAL 
variability. 



5.2 Changing Covering Fraction 

While most of the evidence presented so far favors ioniza- 
tion changes, previous BAL variability studies, including 
Paper 1, have fa vored changing covering fractions over ion- 



ization changes IlLundgren et al.l 20071: IGibson et al. 



lHamann etal] 2008; Krongold et al. 2 O10l : lHall et all 



2008; 



20111 ). 



To investigate this possibility, we consider a simple scenario 
with clouds that have constant ionization and column den- 
sity, but are moving across our line-of-sight. If C IV and Si IV 
have the same covering fraction, then C IV should be just as 
likely to vary as Si IV, which is inconsistent with the results 
of Section 3.1. Further, the change in strength in the two 
lines should be the same, which is con tradicted by our re- 
sults in Fig. [51 (also. IGibson et alJl20TOf ). Hence, this simple 
scenario clearly does not match the results of this and pre- 
vious work. 

A more realistic scenario involves clouds that can 
have different co v ering f ractions in C i v an d Siiv 
llBarlow fc Sargentl Il997l: lHamann et all Il997l. l200ll: 



Ganguly et al. 1999 | : Gabel et al 



2005, 



2006; lArav et all I2008T ). lHamann et al 



Gabel et al.l 

20011 ) 



and 



Hamann fc Sabral ([2004) discuss simple schematics of 



inhomogeneous clouds that could lead to differe nt covering 
fractions in different ions (see fig. 6 in lHamann et al.ll200ll 
and fig. 2 in lHamann fc Sabral 120041 ) . Stronger transitions 
in more abundant ions can have a larger optical depth over 
a larger area in these schematic models. Thus, Siiv may 
trace a different area of the outflowing gas clouds than C IV. 

If t he C iv and Si iv lin es are saturated, e.g. like the 
BALs in lHamann et al.ll2008l . then the strengths of the lines 
would be governed by the covering fractions in those lines. In 
this case, a smaller covering fraction in Siiv would be con- 
sistent with the results of Fig. [2] which shows that Si iv lines 
are generally weaker than C iv. If Si iv is tracing a smaller 
area of the gas cloud than Civ and this cloud is moving 
across our line-of-sight, then Siiv absorption would gener- 
ally be more variable. Furthermore, if the covering fractions 
are different for each ion, then the change in covering frac- 
tion, as well as the fractional change in strength of the ab- 
sorption lines, for each ion can also differ. This is consistent 
with Figs. [gland El 

In Sections 13.3.11 13.3.21 and 13.3.111 we report on three 
cases that show evidence of one BAL, or one portion of a 
BAL, strengthening while another weakens within the same 
quasar. This can be readily explained in a moving cloud 
scenario, for it is possible for clouds at different velocities to 
enter/leave our line-of-sight at different times. If one cloud 
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enters our line-of-sight, while another is exiting, we would 
see one BAL strengthening while another weakens. 



5.3 Conclusions 

The higher variability fractions in Si IV versus C IV B ALs 
and coordinated variabilities between absorption regions at 
different velocities in individual quasars supports the sce- 
nario of global changes in the ionization of the outflowing gas 
causing the observed BAL variability. Furthermore, velocity- 
dependent covering fractions and optical depths could ex- 
plain why in many cases we see variability in just portions 
of BAL troughs, rather than entire troughs varying. On the 
other hand, variability in portions of BAL troughs fits nat- 
urally in a scenario where movements of individual clouds, 
or substructures in the flow, are causing changes in covering 
fractions in the absorption lines. This scenario is also con- 
sistent with the main results of Section 3.1, assuming that 
Si IV has a smaller covering fraction than C IV. 

In reality, changes in ionization and covering fractions 
could both be contributing to BAL variability. In our sam- 
ple, there are quasars in which we observed no variability; 
there are quasars that varied in only one or two narrow ve- 
locity intervals; and, there are yet other quasars with varia- 
tions over a wide range in velocities. It is unlikely that one 
scenario is governing the changes in all of these quasars. Per- 
haps in quasars where the lines are more saturated, the lines 
are not susceptible to small changes in ionization, but can 
easily vary due to covering fraction changes. In other cases, 
where the lines have lower optical depth, a change in ioniza- 
tion can cause large changes over a wide range in velocities, 
possibly masking variations due to changes in covering frac- 
tion. 

There are still some unanswered questions that our re- 
sults from Paper 1, and the current work, raise. In particu- 
lar, in Paper 1, we find correlations between incidence of C IV 
variability and both velocity and absorption strength. How- 
ever, velocity and absorption strength are also correlated. If 
the trend is really with absorption strength, indicating that 
weaker lines, which are less likely to be saturated, are more 
variable, then this favors ionization changes. If the trend 
is with velocity, then the implications are more ambiguous, 
but it would be more consistent with the crossing cloud sce- 
nario. Clouds with higher outflow velocities are more likely 
to have greater transverse velocities as well. 

There are also documented cases of BALs emerging 
where there had hitherto be en no absorption (|Hamann et al.l 
120081 : iKrongold et al.ll2010l ). and in this work, we report on 
a quasar (1303+3048; Section \3. 3. 8[) where a mini-BAL be- 
came a BAL. These scenarios speak to the general com- 
plexity of BAL variability. Previous studies have hypoth- 
esized that different outflow lines may indicate different 
inclinations of our lines-of-sight to the quasa rs and that 
at certain inclinations no absorption is seen (|Elvia 120001 ; 
iGangulv et al J 120011 ). The connection between BALs, mini- 
BALs, and NALs is still unclear. If BAL and mini-BAL 
outflows occur at different inclinations, then perhaps our 
line-of-sight to 1303+3048 goes through a region of overlap 
between the mini-BAL and BAL inclinations. One might 
expect this putative border region between the BAL and 
mini-BAL parts of the flow to be the most turbulent or un- 



stable, and thus the most prone to showing line variability 
caused by structural changes/motions in the flow. 

This work is only the second paper in a 3-part series 
on BAL variability. The next paper will include 1) a more 
thorough exploration of the variability time-scales, with new 
data added to give extensive coverage across week to month 
intervals; and 2) a more complete discussion of the implica- 
tions of variability, e.g., in terms of the size, location, and 
stability of outflow structures. 
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